Draft version December 27, 2012 

Preprint typeset using LAT^X style cmulatcapj v. 5/2/11 



(N 

o 

(N 
O 

Q 

(N 
(N 

< 
Ok 

6 

CO 

o3 



in 

(N 



13 



THEORETICAL EXPLANATION OF THE COSMIC RAY PERPENDICULAR DIFFUSION COEFFICIENT IN 

THE NEARBY STARBURST GALAXY NGC 253 

K. Buffie 1 , V. Heesen 2 , and A. Shalchi 1 

department of Physics and Astronomy, University of Manitoba, Winnipeg, Manitoba R3T 2N2, Canada, shalchi@physics.umanitoba.ca 

and 

2 School for Physics and Astronomy, University of Southampton, Southampton S017 1BJ, UK, v.heesen@soton.ac.uk 

Draft version December 27, 2012 

ABSTRACT 

Diffusion coefficients are usually used to describe the propagation of Cosmic Rays through the 
Universe. Whereas such transport parameters can be obtained from experiments in the Solar System, 
it is difficult to determine diffusion coefficients in the Milky Way or in external galaxies. Recently 
a value for the perpendicular diffusion coefficient in the nearby starburst halaxy NGC 253 has been 
proposed. In the present paper we reproduce this value theoretically by using an advanced analytical 
theory for perpendicular diffusion. 

Subject headings: cosmic rays - galaxies: magnetic fields - galaxies: individual (NGC 253) - galaxies: 
starburst 



1. INTRODUCTION 

Radio continuum observations of the nearby starburst 
galaxy NGC 253 can be used to measure th e distribution 
and transpo rt of Cosmic Ray electrons (see lHeesen et al.l 
I2009all2"bllh . The current understanding is that Cosmic 
Ray electrons are accelerated in t he disk by shock wave s 
induced by supernova explosions ([Reynolds et al.l 12012) . 
The electrons are then transported away over their life- 
time by either convection in a galactic wind or diffusion, 
where the two processes can be distinguished from their 
different typical transport length as a function of elec- 
tron energy. The magnetic field in galaxies is normally 
dominated by the turbulent component, but averaged 
over spatial sizes larger than 1 kpcQ (a typical resolution 
achieved with radio interferometers in external galaxies) 
ordered magnetic field components are detected by their 
linear polarisation from synchroton emission of Cosmic 
Ray electrons spiralling in the magnetic field. Observa- 
tions suggest that galaxies in their disk plane are domi- 
nated by a magnetic field that is parallel to the disk and 
has a spiral pat tern with a constant pitch angle of a few 
tens of degrees (lBeckll2012A . Further away from the disk 
plane at heights larger than 1 kpc the field becomes more 
vertical and the field lines o pen up further into the halo 
(jHaverkorn fc Heesedl2012l ). 

The motion of Cosmic Rays through the interplanetary 
or interstellar space is complicated due to their inter- 
action with turbulent magnetic fields. Such fields lead 
to spatial diffusion of the energetic particles. In addi- 
tion to the turbulent magnetic fields one can also ob- 
serve a non-vanishing mean magnetic field which is in 
the case of the interstellar medium the Galactic mag- 
netic field. This mean field breaks the symmetry of the 
physical system and one has to distinguish between dif- 
fusion along and across the Galactic magnetic field. The 
former process is also known as parallel diffusion whereas 
the latter process is called perpendicular transport. The 
theoretical investigation of Cosmic Ray transport has a 

1 In the present article we measure distances in parsecs (pc). 
Approximately we have lpc PS 3 ■ 10 16 m S3 3.31y. 



long history (see, e.g., Schlickciscr 2002). Whereas the 
propagation of cosmic particles through the solar sys- 
tem seems to be well underst ood (see, e.g.. iShalchi et afl 
120061 : iTautz fc Shalchil l2012h . it is unclear whether ad- 
vanced transport theories can describe particle propaga- 
tion with high accuracy in other systems like the Milky 
Way or external galaxies. At least parallel diffusion 
seems to be a process which can be described very well 
by more advanced diffusion theories. The measured de- 
crease of the abundance ratio of secondary to primary 
Cosmic Ray nuclei as B/C and N/O at kinetic energies 
above 1 GeV/nucleon, implies a variation of the total 
column density as a function of rigidity R as (S wordy 
et al. 1990) An oc R 6 where we have used the paral- 
lel mean free path Ay . This behavior was explained by 
IShalchi fc~ Schlickciscr (2005) who used a nonlinear dif- 
fusion theory. The behavior of low energy Cosmic Rays 
was also ex plained by using an exten sion of quasilinear 
theory (see IShalchi fc BuschinalfeOlOf ). 

More complicated, however, is the process of perpen- 
dicular di ffusion. In the recent years progress has been 
achieved. iShalchil (|2010[ ) developed an advanced theory 
for perpendicular diffusion which contains all known lim- 
its of diffusion theory. This theory is called the Uni- 
fied Non-Linear Transport (UNLT) theory and provides 
the correct subdiffusive behavior for one-dimensional tur- 
bulence in agreement with the so-called theorem on re- 
duced dimensionality (see Jokipii et al. 1993 an d Jones et 
al. 19 98) and computer simulations (see, e.g.. lOin et al.l 
2002). Furthermore, the correct field line random walk 
limit can be derived in agree ment with the n o nlinea r 
field line transport theory of iMatthaeus et al.l (|1995f ). 
The theory discussed here was also tested numerically 
confirming its yalidity for solar wind parameters (see 
ITautz fc Shakai2lmh . 

In analytical theories for energetic particle transport, 
the model of interstellar turbulence enters the corre- 
sponding nonlinear inte gral equation for the perpendicu- 
lar diffusion coefficient. Shalchi et al.l (|2010f) have there- 
fore combined the UNLT the ory with the turb ulence 
correlation tensor proposed by iCho et ail ([2002D based 
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on G oldreich-Sridhar sc aling (see IGoldreich fe Sridharl 
I1995T) . Later, however, iHeesen et al.l ()201lD compared 
the analytical result obtained in the aforementioned pa- 
per with observations of NGC 253 and noticed a major 
difference. Using a multi-wavelength picture of the star- 
burst nucleus they established that the nuclear outflow of 
hot X-ray emitting gas can be collimated by the magnetic 
field lines in the walls that surround the conical outflow. 
The magnetic field lines are aligned with the wall ori- 
entation as would be expected from compression by the 
hot gas in the outflow cone. An annotated figure of the 
multi- wavelength picture can be found in FigureHJ which 
shows the approximate location of the magnetic field fil- 
aments that are the walls of the outflow cone seen in 
projection. We have overlaid contours of the A3 cm 
radio continuum emission at a resolution of 150 pc 
that show the extent of the filaments extending 
along the northwestern outflow cone. The fila- 
ments are also observed at AA 20 and 6 cm at the 
same spatial resolution, which offers to study the 
perpendicular diffusion coefficient as described in 
Section [H 

These observations led [Heesen et al.l ()2011h to 
the finding that the observed perpendicular dif- 
fusion coefficient is a factor of 10 too high in 
comparison to the predicted value by the theory. 

In the present paper we, therefore, revisit the problem 
of perpendicular transport of energetic particles in the 
nearby st arbust galaxy NGC 253. We combine the ap- 
proach of IShalchi etaLl (|2010f ) with more realistic parti- 
cle and turbulence parameters to compute the perpen- 
dicular diffusion coefficient more accurately. 




Figure 1. Three-colour composite, multi-wavelength view of the 
central region in NGC 253. Red, green, and blue indicate Ha from 
IWestm oaucttc at al. (2011), A20 cm radio continuum, and Chan- 
dra soft X-ray (M. Hardcastle, priv. com.), respectively. Contours 
show A3 cm radio continuum at 150 pc resolution with the first 
contour at 3 X the r.m.s noise level and each consecutive contour 
increasing roughly by a factor of two. An id entical figure withou t 
the annotations and contours can be found in Heesen et al. (2011). 



a review). Very recently an advanced theo ry for perpen- 
dicular transport has been developed (see lShalchll2010l) 
which shows good agreement with computer simulations 
performed for the pa rameters of the solar wind (see 
iTautz fc Shalchil 120111 ). This theory is called the Uni- 
fied Non-Linear Transport (UNLT) theory and provides 
a nonlinear integral equation for the perpendicular dif- 
fusion coefficient k±. According to this equation k± de- 
pends on magnetic turbulence described by the so-called 
magnetic correlation tensor. Therefore, the understand- 
ing of turbulence is crucial for the analytical description 
of cross field diffusion. A model for the magnetic cor- 
rel ation tenso r in th e interstellar medium was proposed 
by iCho et aT] (120021) which is based on the Goldreich- 
Sridhar model (IGoldreich fc Sridharlll 995). For this spe- 
cific turbulence model, IShalchi et al .1 (|2010D derived the 
following equation 

2.?V + 1 

X x 2 y 2 + l x 2 y 4 X ll /X± + y 2 + 4A,, A ± /(3/) 2 1 ' 

where we have used the perpendicular mean free path 
which is related to the perpendicular diffusion coefficient 
via X± — 3kj^/3. According to Eq. ([1]) the parameter 
Aj_ depends on the parallel mean free path A|| = 3k\\/v 
(v is the particle speed) and, therefore, there is a mutual 
influence of parallel and perpendicular transport. The 
parameter a 2 in Eq. ([1]) is an order one constant related 
to the probability that the particle is scattered away f rom 
the magnetic field lines (see, e.g.. lDosch et al.|[2009D . In 
the present paper we set a 2 = 1 for simplicity. The per- 
pendicular diffusion coefficient is also controlled by the 
characteristic length scale of the turbulence I, and the 
ratio Eb = 5B/Bq. The latter ratio is a measure for 
the turbulent magnetic field strength SB with respect 
to the mean field Bo. The mean field is Bq defined as 
the spatially averaged magnetic field vector on a scale 
corresponding to the length I, which is usually related 
to the correlation length of magnetic turbulence. The 
latter parameter is a characteristic scale for the correla- 
tion of magnetic fields at two different positions in space. 
Often this scale is directly proportional to the so-called 
bendover scale which denotes the turnover from the in- 
ertial range of the turbulence to the energy range. Fur- 
thermore, this scale is smaller than the largest scales on 
w hich turbulence can be observed. 

IShalchi et al.l (|2010f ) have solved Eq. (fT]) numerically to 
obtain values for the perpen dicular diffusion coe fficient 
in the interstellar medium. IHeesen et al.l ()2011D . how- 
ever, compared the theoretical results with observations 
and found a disagreement. Below we, therefore, perform 
a parameter study to compute the perpendicular diffu- 
sion coefficient and compare it with observations. We 
will show that for the correct choice of input parameters, 
Eq. ([1} is indeed able to reproduce the aforementioned 
observations. 



2. THEORETICAL DESCRIPTION OF 
PERPENDICULAR DIFFUSION 



The analytical description of diffusio n across the m ean 
magnetic field has a long history (see iShalchil 120091 for 



3. OBSERVATIONS 



The different parameters entering Eq. (p} as well is the 
perpendicular diffusion coefficient itself can be obtained 
from radio continuum observations of the nearby galaxy 
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NGC 253 as described in lHeesen et~aH (|2009at 12011ft . In 
this section, we briefly summarize the techniques they 
have been using to quantify the Cosmic Ray transport. 

Radio continuum studies can be used to measure the 
length scale of Cosmic Ray diffusion. The time scale of 
Cosmic Ray electrons is determined by loss processes by 
which the electrons are losing their energy, such as syn- 
chrotron and inverse Compton radiation. The diffusion 
coefficient can be then calculated by 

Hobs = £diff/ T J (2) 

where Ldiff is the diffusion length scale and r is the elec- 
tron lifetime. Depending on the magnetic field structure, 
the observed diffusion coefficient is either along (fty) the 
magnetic field or perpendicular to it (k±). 

Ideally, one can measure the distance to the star- 
formation sites, where the Cosmic Rays are accelerated 
and injected into the interstellar medium, to obtain the 
transport length scale. This is for instance the case when 
observing galaxies in so-called edge-on geometry, where 
the observer is located in or close to the disk plane of 
the galaxy (the inclination angle is close to 90°). Cos- 
mic Ray acceleration in supernova remnants is confined 
to the relatively thin disc plane, where the formation of 
massive stars happens. The geometry in this case is thus 
simple: Cosmic Rays are transported away from the star 
formation sites over their lifetime and their transport 
length scale is equal to the vertical electron scaleheight. 
Typical scaleheights o f galaxies are 1.8 kpc at observing 
wavelengths of A6 cm (|Krausell2009|) , where longer wave- 
lengths have larger scaleheights and shorter wavelengths 
smaller ones. Typical Cosmic Ray electron lifetimes in 
the interstellar medium are between 1 and 10 Myr. This 
basically determines the typical parallel diffusion coef- 
ficient in the ISM, where it is assumed that the halo 
magnetic field is opening up from the disk parallel to a 
vertical direction with increasing distance from the disk. 
Therefore, the vertical diffusion is predominantly along 
the magnetic field lines hence allowing us to measure K|| . 
IHeesen et ahl (|2009bl ) confirmed this for NGC 253 finding 
the characteristic X-shaped halo magnetic field with sig- 
nificant vertical compon ents that is observ ed in a number 
of nearby galaxies (e.g.. iSoida et al.|[2011|) . 

To measure the perpendicular diffusion coefficient is 
more difficult, because it requires to have variations of 
the Cosmic Ray distribution perpendicular to the mag- 
netic field orientation. The observed radio continuum 
emission is mostly smooth in the disk-halo interface and 
filamentary structures are rare. Numerical MHD simula- 
tions suggest that the disk ha lo interface is dominated b y 
filamentary magnetic fields (jBreitschwerdt et al.ll2012n . 
but line-of-sight confusion and limited spatial resolution 
hampers their detection. An exception are starburst 
galaxies, where the spatially concentrated star forma- 
tion activity results in exceptionally high radio contin- 
uum surface brightness, allowing us to employ high spa- 
tial resolution (w 100 pc). The spatially concentrated 
star formation can result into outflows of hot X-ray emit- 
ting gas in a galactic wind. The magnetic field is then 
concentrated and amplified by expansion of the hot gas 
until a pressure equili brium is reach e d. T his very spe- 
cific geometry allowed IHeesen et al.l (|2011| ) to measure 
the perpendicular diffusion across the magnetic field in 



the walls of the nuclear outflow cone in NGC 253. 

The latter authors found the width of radio continuum 
filaments weakly dependent on frequency and thus elec- 
tron age, which can be interpreted as perpendicular Cos- 
mic Ray diffusion. We note that they were able to mea- 
sure the orientation of the ordered magnetic field along 
the filaments as it would be expected for a compression 
due to the hot X-ray emitting gas outflowing from the nu- 
clear starburst. From the observations of particles with 
a magnetic rigidity of R = 3 x 10 12 Volt (equivalent to 
an electron energy of 3 GeV), they found a perpendicular 
diffusion coefficient of kj_ = (2.6 ± 0.6) x 10 28 cm 2 s _1 . 
For the parallel diffusion coefficient one can assume 
«ll = 1.0 x 10 29 cm 2 s _1 (again for an electron energy of 
3 GeV) , as measured from the diffusion of Cosmic Ray 
electrons from the disk into the halo along the verti- 
cal halo magnetic field (|Heesen et al.ll200 9a). It is im- 
portant to measure the diffusion coefficients at roughly 
the same electron energy, because they are energy de- 
pendent. We have reduced the measured diffusion co- 
efficient of IHeesen et al.l ()200 9a) to account for a possi- 
ble contribution of convection. The southwestern halo 
is different from the northeastern one as indicated by 
the dependence of the electron scaleheight on the elec- 
tron lifetime and by the different amount of extra-planar 
gas that is more abundant in the northeastern connec- 
tive halo. However, we can not rule out the contribution 
of convection of Cosmic Rays by the disc wind and thus 
have lowered the diffusion coefficient assuming that both 
convection and diffusion contribute equally in the south- 
western halo. 

Equation |lj depends on turbulence properties, namely 
the ratio of magnetic fie l ds SB /Bo and the scale I. Ac- 
cording to IHeesen et al.l ()2011h the equipartion estimate 
for the total field strength is B tot — 46 /iG. The ordered 
field strength -B rd,_L in the sky plane can be computed 
from the observed deg ree of polarization p of the syn- 
chrotron emission (e.g. JBeck fc Krausd 12001 : 

p = p (l + l q 2)/(l + 3q 2 + f q 4 ), (3) 

where po is the intrinsic degree of polarization (po = (3 — 
3a) /(5 — 3a)) and q is the ratio of the isotropic turbulent 
field -Bturb to the ordered field -B rd,± in the plane of the 
sky. For a non-thermal radio spectral index of a = — 1 
and a polarization degree of p — 0.21 we obtain -Bturb = 
41 /iG and i? rd.± = 21 /iG. For this particular case, the 
ordered magnetic field lies in the sky plane so that we find 
for the ordered magnetic field strength Bo — 21 /xG. The 
fluctuations in the magnetic field strength are equal to 
the turbulent field field strength (SB = -Bturb) and thus 
we find for the the ratio SB /Bo w 2. We note that this 
is an upper limit as our spatial resolution of 150 pc that 
was available for the polarization measurements may be 
not be enough to resolve the filamentary magnetic fields. 
High resolution A 20 cm maps show a filament width of 
only 40 pc in radio continuum emission (see Figure [I}, so 
that our measured ordered magnetic field strength B rd is 
only a lower limit. We have thus studied the implication 
of a lower value of SB /Bo — 1 to indicate how it would 
change the theoretical expectation for the perpendicular 
diffusion coefficient. 

More difficult to estimate is the scale I. The width of 
the cone walls in which the magnetic fields are confined 
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is equal or less than 40 pc (jHeesen et al.ll20lH ). This sug- 
gests that the upper value cannot be larger than about 
I ~ 50 — 100 pc. However the value for I is very uncertain 
and, therefore, we compute the perpendicular diffusion 
coefficient for a whole range of correlation lengths. Beck 
(2007) for instance suggested that the largest scales of 
turbulence are in the order of 10 — 100 pc. Such largest 
scales, however, are not necessarily equal to the tur- 
bulence correlation scale. Actually they can be seen 
as the maximum o f the scale I. Sometimes (see, e.g., 
iShalchi et al.l [2009) it is assumed that the correlation 
length of interstellar turbulence is 1 pc but it could also 
be shorter. Therefore we compute the perpendicular dif- 
fusion coefficient for 1 pc < I < 1000 pc to explore the 
values of I which lead to agreement between theory and 
observations. 




10° 10' 10 z 10 3 



I in parsec 



Figure 2. The perpendicular diffusion coefficients of Cosmic 
Rays in the nearby starbust galaxy NGC 253. We show the 
theoretical perpendicular diffusion coefficient for different length 
scales I. For the parallel diffusion coefficient we have used K\i = 

1.0 X 10 29 cm 2 s _1 as proposed by IHeesen et al.l d2009al') . The 
theoretical values were calculated for 5B/Bq = 1 (dashed line) 
and SB /Bo = 2 (solid line). Th e dot represents dif fusion coef- 
ficient from the observations (sec Hccscn ct al. 2011), where we 
have K± = (2.6 ± 0.6) X 10 28 cm 2 s _1 at I = 50 pc. 



4. REPRODUCING THE OBSERVATIONS 

In the following we solve Eq. (JT|) numerically for differ- 
ent turbulence scales I. To relate the mean free paths to 
the spatial diffusion coefficients we employ A|| = 3k\\/v 
and Aj_ = 3k±/v. In the latter relations we replace the 
particle speed v by the speed of light c since we deal 
with relativistic particles. As described in the previous 
sections we set k« — 1.0 x 10 29 cm 2 s _1 and a 2 = 1 in 
our equation for the perpendicular diffusion coefficient. 
We have calculated the latter transport parameter for 
SB /Bo = 1 and SB /Bo = 2. Our theoretical results are 
shown together with the observations in Fig. [2] 

The observational value of the perpendicular diffusion 
coefficient is (2.6 ± 0.6) x 10 28 cm 2 s~ 1 for a correlation 
length of approximately 50 pc. The theoretical perpen- 
dicular diffusion coefficient is calculated for different val- 
ues of the correlation length. For a correlation length of 
50 pc we found a K^which is very close to the observa- 
tions depending of the value of 5B/B - We note that 



the perpendicular diffusion coefficient is only weakly de- 
pendent on the turbulent correlation scale above a scale 
of approximately 10 pc. The observed diffusion coeffi- 
cient seems to rule out correlation lengths smaller than 
lpc, but does not constrain very well the upper limit. 
It is reasonable to assume that the correlation length is 
smaller than the largest sca les which are in the order of 
10 - 100 pc (see lBeck|[2007h . We note that the observed 
value of the perpendicular diffusion coefficient is quite 
well reproduced by the measured value for SB/B$, which 
is in the range between 1 and 2 as argued in Section |3J 

5. SUMMARY AND CONCLUSION 

A fundamental process in astrophysics is the propaga- 
tion of Cosmic Rays through the Universe. Whereas it 
seems that we understand the motion of e nergetic parti- 
cles in the solar system very well (see, e.g.. IShalchi et al.l 
l2006HTautz &: Shalchi| [2012). it is still not clear how the 
transport parameters look like in the interste l lar sp ace of 
our own or external galaxies. Sh alchi et al.l ([2010D have 
calculated the perpendicular diffu sion coefficien t for in- 
terste llar turbulence parameters. IHeesen et al.l (2009a, 
1 2 llf) obtained measurements for the parallel and per- 
pendicular diffusion coefficient in the nearby starbust 
galaxy NGC 253. They comp ared thei r resul ts with the 
theoretical values presented in IShalchi (|2010f) and found 
one order of magnitude between the two results. How- 
ever, the perpendic ular diffusion coefficient obtained in 
IShalchi et ail (l2010l) was obtained for a different set of 
parameters. Especially the parallel diffusion coefficient 
used in the latter paper does not agree with the value ob- 
tained by the observers. Thus, in the current paper we 
combined the approach of IShalchi et al.l (|2010f ) with pa- 
rameter sets tailored to the observations of NGC 253, 
such as the parallel diffusion coefficient and magnetic 
fields strengths. Furthermore, we computed the param- 
eter k± for different values of the correlation length I 
as this value is known only to a small degree of cer- 
tainty. Our results now largely reconcile the observed 
values with the theoretical predicted ones as visualized 
in Fig. [2] For correlation lengths between 5 and 100 pc, 
the theoretical value and the observed one agree within a 
factor of 2. This agreement is quite good given the diffi- 
culties in the observational measurements of diffusion co- 
efficients for Cosmic Ray electrons from radio continuum 
observations in galaxies. We can confidentially exclude 
correlation lengths smaller than 1 pc, which would result 
in too small perpendicular diffusion coefficients. 

We have done this study for only one galaxy, because 
so far measurements particularly of the vertical diffu- 
sion coefficient are very sparse in galaxies. However, the 
agreement of the observations with the theoretical pre- 
diction lends to some degree support both to the theoret- 
ical description of Cosmic Ray diffusion as described in 
this paper and to the observational attempts to measure 
it from the Cosmic Ray electron distribution. 
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